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ABSTRACT 


Work in the area of mechano-caloric phenomena has been initiated during 1988 with 
startup in the Summer 1988 period. The ideal system utilizing He II "super-phenomena" is 
modeled readily, within the frame of thermodynamics energetics, using the concept of an ideal 
superleak. The real system however uses porous media of non-ideal pore-grain ingredients. The 
early phase of experimental and related modeling studies is outlined for the time period from 
Summer 1988 to the end of 1988. 


1 



TABLE OF CONTENTS 


Page 

1. ABSTRACT 1 

2. INTRODUCTION 3 

3. 1988 PROGRESS AND DISCUSSION 4 

4. REFERENCES 17 

APPENDICES 

A. ONSET OF VAPORIZATION ASSOCIATED WITH COUNTERFLOW 

IN POROUS MEDIA 18 

B. CRITICAL TRANSPORT PARAMETERS FOR POROUS MEDIA 

SUBJECTED TO COUNTERFLOW „ . 

24 

C. CHARACTERISTIC PARAMETERS RELATED TO SUPERCONDUCTOR- 

COOLANT INTERACTION , n 

40 

D. ABSTRACT OF M.S. THESIS OF P. KHANDHAR 59 


2 



2 . 


INTRODUCTION 


The project has been started up with the goal of obtaining realistic performance values of 
a mechano-caloric component The mechano-caloric effect has been documented in the early 
literature on the subject 1 . The effect is one particular phenomenon out of several peculiar pro- 
perties of "superfluids". The spectacular property of "superfluidity" has been reported in 1938 by 
Kapitza and Allen and Misener. The latter authors used a bundle of fine stainless steel wires, 
packing it into a thick walled, low-conductivity alloy tube ("German silver"), and drawing the 
system down to a small diameter. The result had been a set of channels between the wires with 
not more than one micro-meter in diameter. Similarly, Kapitza experimented with annuli 
arranged to have a distance, between walls, on the order of magnitude 1 micro-meter or less. 

The velocity of flow of liquid He II (i.e. Helium-4 below its lambda point), turned out to be 
nearly independent of an externally applied pressure head. This unusual property has been 
referred to as "superfluidity". If flow rates are unaffected by the driving force, the flow resis- 
tance commonly observed does not exist. Hence the term superfluid liquid Helium appears to be 
justified. 

The next most important "new effect" found with liquid He II has been the ther- 
momechanical effect: a difference temperature (AT), applied externally by suitable means, 
causes a pressure difference (AP). The first observations have been reported by Allen and 
Jones 2 . The thermomechanical pressure difference in these early experiments has been con- 
verted to kinetic energy associated with a "fountain" of liquid He II. Hence the term fountain 
effect has become popular as well. The presently accepted two-fluid model macroscopics postu- 
lates a simple description involving two separate fluid components (normal and superfluid). The 
additional postulate of "zero chemical potential difference" (d|i = 0) for this system leads to the 
famous H. London equation (e.g. Wilks, 1967) 3 . The result of thermostatic conditions is 
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(3P/3T)^ = pS 


( 2 . 1 ) 


(p liquid density, S entropy per unit/mass of liquid). 

As noted by Wilks (op. cit. above), there is a corollary to the thermomechanical experi- 
ments: an application of an externally applied pressure difference (head) may lead to heating or 
to cooling. Equation (2.1) is rewritten in the form 

(3T/3P) = l/(pS) 

P (2.2) 

When the pressure increment is dP > 0, the incremental temperature change has to be dT > 0, i.e. 

heating is observed. When the pressure drops however, (dP < 0), there is a decrease dT < 0. 

This implies mechano-caloric cooling. This effect was observed first by Daunt and Men- 
delssohn 4 . 

The two-fluid interpretation is as follows: only superfluid, devoid of viscosity, will flow 
through the narrow fluid passages cited above. The consequence is an enrichment of normal 
fluid on the other side of the passages. An increase in normal fluid "concentration" means an 
enhanced temperature. If the sign changes such that superfluid flows into the fluid space con- 
sidered, a drop in T occurs. Thus, cooling may be accomplished in a refrigerator, e.g. the He II 
"vortex fridge" (Staas-Severijns,1969) 5 . 

3. 1988 PROGRESS AND DISCUSSION 

The preceding section points out the potential use of the He II "supereffects" in unique 
refrigeration devices called He II vortex refrigerators. Their components are the pump (possibly 
a thermomechanical effect device), and a mechano-caloric component. During the work done up 
to this time, it became clear to us that a careful plug selection is possible only when irreversibili- 
ties along with most favorable thermodynamic parameters are known. The following topics are 
considered: phenomenology of porous media, metal porous plugs at zero net mass flow 
(ZNMF), and ceramic plugs at ZNMF. 
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Phenomenology of porous media . Figure la is a schematic diagram of the ideal and real pres- 
sure difference (static case) of a He II thermomechanical pump (fountain effect pump, FEP). 
Figure lb schematically presents a similar median o-caloric component (MCC) diagram. Figure 
la indicates that the ideal superleak case (ISL) is an upper bound for an extended temperature 
difference. Irreversibilities at large AT cause a saturation in the static AP. 

Figure lb indicates that the ISL temperature difference predicted is an upper bound to the 
real AT when irreversibilities cause a less efficient change of state than the ISL case. Figure 2 is 
a plot of the pressure increase of the ideal pump (ISL) for various initial bath temperatures. Fig- 
ure 3 is a plot of the ideal lowering of the temperature in a mechano-caloric component (ISL 
case). The consequence of the calculations and of early data is a need to get experimental 
evaluations of particular plug materials for each case. Furthermore, it is necessary to obtain a 
unique measure of the liquid throughput capability in a real plug. The quantity chosen is the 
Darcy permeability of a porous medium. The so-called Zero Net Mass Flow (ZNMF) technique 
has been used in order to characterize the plugs employed in the research. Details of these stu- 
dies are outlined in the Appendix sections A and B. 

Metal porous plugs at ZNMF . The present metal systems are in a permeability range below the 
more thoroughly documented plug range of vapor-liquid phase separators. The latter have been 
described in earlier work, e.g. the Ph.D. thesis of Dr. Sidney Yuan (UCLA, 1985). Figure 4 is a 
characteristic length diagram from the M.S. thesis of P. Khandhar. The abstract of the thesis is 
attached as Appendix D. The present work suggests that metal plug manufacturers have 
improved their fabrication capabilities. They produce well-defined plugs in the Darcy permea- 
bility (order of magnitude) range of k - 10 -10 cm 2 , while several years ago KH-IO -9 cm 2 was 
available only. 
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Fig. 1. Schematic diagram of (absolute) temperature difference (AT) and pressure difference 
(AP) behavior. 

a. fountain effect pump (FEP) 

b. mechano-caloric component (MCC) 
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Figure 4 contains packed bed ceramic plugs of Singsaas (Ref. 6). These plugs have been 
considered less desirable for space conditions than rigid plugs whose particles are set. In addi- 
tion some mechanical rigidity without excessive exposure to forces appears to be useful. While 
the metal plugs are documented in the Appendix sections A and B, the ceramic experiments 
have been in a rather early stage. 

Ceramic plug at ZNMF . We have used two ceramic plug systems. One has been cooled down 
only to the liquid nitrogen range. The other one has been subjected to liquid Helium conditions. 

Both plugs are units manufactured by Coors Ceramic Corp., Golden, Colorado. The 
ceramic components designed for filtration purposes have the shape of a cylindrical cup, open at 
one end. The geometry does not have a perfect cylindrical shape because of firing constraints 
after initial shaping. Two cups are joined together as shown schematically in Figure 5. 

A phenolic ring (PR) has a dual purpose. First, for low heat leakage it has a small non- 
porous area (compared to the total cup area available for porous media flow); second, the various 
instrument wires are routed through the ring. Instrumentations include an internal heater unit 
(IHU), a pressure transducer (PD) (Siemens KPY type), and a carbon thermometer (C-TH). 

Pressure transducer characteristics have been determined*. An example is Figure 6. 
Various results obtained appear to be consistent with early manufacturer data 7 . 

Results obtained with the second sup set are presented in Figure 7 as pressure transducer 
response as a function of the heat input during the ZNMF runs. These data are for rather high 
heat inputs (compared to the critical rates, estimated for these rather fine plugs). When the plug 
unit had been warmed up, after the He II runs, it turned out that there were cracks visible after a 
thorough room temperature inspection. Areas detected have been covered with epoxy. After a 


*Dr. A. Eisner has dynamically contributed to the early pressure transducer studies, and his effort is acknowledged 
with appreciation. 
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Fig. 5. Ceramic porous filter cup geometry, schematically. PR phenolic ring, IHU internal 

heater unit, C-TH carbon thermometer, PD pressure transducer. 
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second set of runs, the same type of cracks had been formed again. The room temperature 
inspection showed additional local domains where ceramic material had fractured on the cup 
surface. Chips, found below the plug/cup system, have been an indication of the brittle nature of 
the ceramics and related lack of mechanical reliability. It is noted that this type of potential frac- 
ture danger is inherent in the recently discovered metal oxide superconductors with high transi- 
tion temperature (HTSC). The assessment of HTSC has been the topic of a recent seminar on 
heat transfer related to stability of superconducting magnets. This work is outlined in Appendix 
C. 


Plug data have been collected for various heat inputs Q into the porous plug chamber at 
different bath temperature. Examples of data are presented in Figs. 8 and 9*,(A total plug cross 
section). Q/A versus temperature difference is presented in Fig. 8. The data indicate a linear 
region in the ZNMF mode at low Q. Figure 9 represents the data in the AT range from 0 to 100 
mK. The data appear to be mostly in the non-linear range. 

Acknowledgments. The contributions of Dr. A. Eisner and of our collaborators David Y. Ono 
and F. Afifi are acknowledged with thanks. 
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Figure 8. Heat 
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APPENDIX A 

ONSET OF VAPORIZATION ASSOCIATED WITH COUNTERFLOW IN POROUS MEDIA 


T. H. K. Frederking, P. K. Khandhar, P. Abbassi 
University of California, Los Angeles, CA 90024 


Porous sintered metal plugs have been subjected to counterflow. A step input in 
power is applied to a heater inside a He I chamber which is connected via the 
plug to an external He n bath. Phenomena of phase transitions are analyzed on 
the basis of thermograms obtained (chamber liquid temperature vs. time data for 
relatively large diameter duct$ there is He II superheating across the first order 
transition phase boundary, of the vapor pressure curve. For high heat inputs, 
extended vapor domains are produced in the chamber for temperatures above 
the homogeneous nucleation limit. 

INTRODUCTION 

He II fluid management and thermal control involving porous media devices have been receiving increased atten- 
tion. Furthermore, the use of narrow coolant channels in magnets produces a special type of "porous" structure. 
The operation of these systems in metastable liquid states, beyond the equilibrium state phase boundaries has been 
the topic of several recent studies [1-6]. The present studies have the purpose of extending the geometries covered 
to include porous media. The following phenomena are encountered: 

Various types of superheating of He II and He I; 

Phase transition involving metastable liquid below the saturated liquid lambda point pressure; 

Phase transitions leading to high temperature excursions into the vapor domain. 

PHASE TRANSITION PHENOMENA INVOLVING FLUIDS (He 4 ) 

The studies cited in references 1-6 have clarified details of phase changes during heat flow in duct systems. The 
findings are briefly outlined, emphasizing pressures (P) below P^m, where P^t is the lambda point pressure of 
saturated He II at liquid He n-He I-vapor equilibrium. For P > Px^*,, there has been no superheating of He II to 
metastable He n states beyond the lambda line during heating. Instead He I domains are formed when the liquid is 
heated to T > T lambda (P). At P < Px*, t , the following phenomena have been observed [1-5]: During a finite dura- 
tion of time (t). He II may be heated up, in a local "hot” domain across the first order phase transition taking place at 
the vapor pressure curve. This involves the latent heat of vaporization on the order of 20 J/(g K). 

According to theoretical and experimental results, the metastable He II cannot be superheated (in a He II state) 
beyond the lambda line extension from P^m, to low P. This line is very steep in the P(T) phase diagram, and T 
lambda is near 2.172 K. 

Metastable He I may be formed at "high T" states, beyond the extended lambda line, in heated domains. From this 
metastable He I, vapor bubbles of nucleate boiling may be formed ("triple-phase nucleate boiling", ref. 7). The 
nucleate boiling mode is quite efficient in keeping the temperature at a heater at low T, close to the Kapitza- 
resistance-dominated solid temperature. 

Metastable He I, upon heating at a high heat flux density q, may be converted into "hot" vapor domains at which 
point the homogeneous nucleation line is traversed during the temperature excursion T(t). The process leads to low 
film boiling rates ("triple-phase film boiling") of heat removal, or high local temperatures at the heater. 

*) Presented as paper C4.5 at the 1988 Intern. Cryog. Eng. Conf. Southampton, 
July 12-15, 1988 . 
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EXPERIMENTS 


The zero net mass flow (ZNMF) mode of the two-fluid model is established. Heat flow takes place during 
counterflow of normal fluid and superfluid in the porous plug. One side of the plug is connected to a heating 
chamber, the other side to an external He II bath. A transient technique is employed, imposing a step input in power 
externally. The resulting thermograms T(t), e.g. Figures 1 and 3, are considered subsequently. 

The porous media (inset of Figure 2) are sintered stainless steel plugs with axis oriented vertically. The chamber 
walls are made out of phenolic (Fig. 2, outer diameter 1.5 cm, overall length 3.6 cm). The heater is wound bifilarly 
to a resistance near 100 ohms using constantan wire (75 pm diameter). The plugs have the following dimensions: 
PKS: diameter D * 6.4 mm; thickness L = 0.75 mm; PSM 05: d = 7.9 mm; L = 1.75 mm; PSM triple: D = 3.23 
mm; thickness of the single plug: L = 1.6 mm, resulting in a total thickness of 4.8 mm. 

THERMOGRAMS AND POROUS PLUG PHASE TRANSITION RECORDS 

The T(t) curves observed show single phase regimes, phase change onsets, and evolution of domains involving 
triple-phase modes of heat flow. 

(a) Single-phase He II modes 

Two principle regimes of transport are observed in the T(t) thermograms. At low step inputs in power, the tempera- 
ture rises monotonically with dT/dt > 0, toward a final steady temperature difference with respect to the outer bath 
(Fig. 1). The normal fluid flows in the laminar regime with a superficial heat flux density, q per total plug cross sec- 
tion A, i.e. Q s qA. The q value is related to a superficial normal fluid velocity (approach velocity) of 

Vno = KnlgradP T |/nn (1) 

(T| n shear viscosity of normal fluid, k„ normal fluid Darcy permeability of the sintered plug). The driving pressure 
difference, established at very small velocities, is London’s thermomechanical difference APt. 

APt = pSAT ; ATcT (2) 

(p liquid He n density). The single-phase He II mode is terminated at the critical value AT C . 

AT C = const(p,/p„)L c 3 p | AP t | /(Lq n 2 ) (3) 

The constant is found to have a particular value for each plug. The characteristic length is L* = (K n ) Vi . 

The non-linear regime is established above AT C . There is a quantitatively different rate of temperature rise (dT/dt) 
initially. At long times, q 0 = Q/A is considerably less than the wide duct value for the same AT. 

(b) Phase transition phenomena at moderate heat input 

The onset of a phase change is marked in the present thermograms by a change of the temperature-time derivative 
from positive to negative values, e.g. Fig. 1 (runs 2 and 3). In principle, such a change may be produced via the 
mechano-caloric effect. However, the ZNMF mode in the relatively wide pore sized range covered appears to lack 
such a behavior with dT < 0. The vaporization phenomenon causing the relative lowering of the liquid chamber 
temperature to dT/dt < 0 appears to utilize the existing excess enthalpy of heated He II in the chamber. This 
implies a superheating of He II across the vapor pressure curve phase diagram (Fig. 6). Related onset values are 
given in Table I. 
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Table I. Thermogram data Qt(U) for onset of vapor formation 
[plug Pk 10-S-02-6.6x0.75; Bath: T = 2.6 K] 

Qt Time t, sec 
0l71 044 

0.153 0.556 

0.123 1.04 

0.115 1.48 

Above Q t , temperature oscillations are observed in a discrete Q-range during the step power input No final steady 
value is found. Instead, the non-linear oscillations persist around a mean T value. These patterns are reminiscent of 
previous observations, however there have been pattern variations from one plug to the next The onset time ti 
appears to change with Q, similar to Q t (tt) in Table I. 

(c) Phenomena at high heat input 

At high q, oscillations may occur, as in (b). However, there is only a restricted duration of liquid T variations at 
temperatures not far from Tx. Figure 3 is a high-q thermogram. After the initial pattern of type (b), the liquid tem- 
perature rises sharply to T above the thermodynamic critical temperature T c = 5.2 K (Fig. 6). 

Figure 4 is a plot of the thermal conductance Q/AT derived from the thermogram (Fig. 3), assuming steady heat 
flow. For the plug area on the order of 0.1 cm 2 , the ratio of qJAT attains the order of 10 W/(cm 2 K). This high 
value indicates that heat flow through the plug is not deteriorated seriously prior to the T-excursion. The rapid drop 
of Q/AT in Fig. 4 signals film boiling conditions. 

Figure 5 shows vaporization onset conditions, corresponding to (b), as Q versus AT. The AT values do reach into 
the metastable He II domain signaling superheating across the vapor pressure curve. The inset of Fig. 5 is a plot of 
AP, recorded by a differential P-transducer on plug (PA/FA PSM 01S-tr-3.2xl.6). The AP-records show an 
insignificant increase in pressure in the non-linear regime (a) of a few millibar. 

The phase diagram P(T) contains the low-P extension of the lambda line and the homogeneous superheat limit. 1 A 
few data points of (b) at and near vaporization onset are included in Fig. 6. Further, several high AT values are 
marked. It is noted that there is no chance to find metastable liquid, for an extended period of time, in the region 
beyond the homogeneous nucleation limit 

CONCLUSIONS 

From the present porous media thermograms it is concluded that our plugs display phase change phenomena con- 
sistent with phenomena observed in larger duct systems connected by the counterflow mechanism to a heated He II 
reservoir. 

Acknowledgment. After initial seed support by the UCLA Senate Research Committee, the investigations have 
received partial support by NASA Ames Research Center. We acknowledge this along with the help and input of 
the Cryogenic Lab. group, in particular, Messrs. F. Afifi, Herbert Simanjuntak and David Ono. 
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Fig. 1. Thermograms far step power input to heater ; (heater current I); V thermometer 


voltage (in arbitrary units A.U.); Plug PK K 10 S 02-6.4 x 0.75 ; 
A. Low power* imposed as step Input in the single-phase He II regime ; 



Fig. 2. Heat flow rate Q versus AT for long tines (t) [steady heat flow]; plug 
PA/FA PSM 01 S-tr3.2 x 1.6 ; 

Inset : Schematic of porous plug system for aero net mass flow ; H heats’ ; 
PH phenolic sleeve ; PP sintered porous plug ; TH carbon resistance thermo- 
meter. 
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CRITICAL TRANSPORT PARAMETERS FOR POROUS MEDIA SUBJECTED TO COUN- 
TERFLOW 

T. H. K. Frederking, F. A. AM, D. Y. Ono 
University of California, Los Angeles, CA 90024 

Experimental and theoretical studies have been conducted to determine critical 
parameters at the onset of non-linear counterflow in He II below the lambda point 
of Helium-4. Critical temperature differences have been measured in porous 
media for zero net mass flow and for Darcy permeabilities in the order of magni- 
tude range from 10 -10 to 10 -8 cm 2 . The normalized critical temperature gra- 
dients, which covered the liquid temperature range of 1.5 K to the lambda tem- 
perature, are found to vary with T proportional to the ratio of the superfluid den- 
sity to the normal fluid density. This liquid temperature dependence appears to be 
consistent with duct data which are limited at low temperature by a Reynolds 
number criterion. 


**) Presented at the 1988 Space Cryogenics Workshop 
Frascati , July 18-19, 1988. 
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Introduction 


There have been numerous reports on "critical" transport parameters at various flow tran- 
sitions of superfluid liquid He II below the lambda point. The two-fluid model, a most success- 
ful concept, has permitted considerable progress in the phenomenological description of the vari- 
ous transport phenomena in He n. According to the basic two-fluid postulate there is superfluid, 
the quantum liquid capable of motion without friction, and normal fluid carrying entropy S and 
heat Q. Both transport rates (S and Q) are subject to dissipation by friction processes. 

One set of critical parameters is concerned with superfluid, in particular with its zero 
resistance flow limits at superfluid critical velocities. Another set of criteria addresses questions 
of flow transitions in heat and entropy convection. Because of viscous forces impeding normal 
fluid motion, there is a finite resistance in all regimes of normal fluid motion. Thus questions 
come up in many application areas concerning the prediction of critical rates. It is the purpose 
of this work to present results of critical temperature difference parameters for porous media 
subjected to counterflow (ZNMF). Despite the restriction to heat flow conditions, it is noted that 
mass effects associated with superfluid (e.g. superfluid density p s ) may enter critical parameters. 
Therefore, a brief discussion of several critical rates is presented first. Subsequently, the discus- 
sion will focus on the experiments. From the data sets collected, we obtain criteria inspected in 
their temperature behavior. Finally, the present experimental data are compared with other 
results and conclusions are presented. 

Critical quantities and major transport modes 

Critical Superfluid Flow. One set of criteria relates to the transition from resistance-less 
"superflow" to resistive flow. The "superflow" cessation has been formulated in terms of quan- 
tum constraints for He II. At the onset of resistance, quantum vortices are generated causing 
vortex shedding into the liquid. The kinematics of these Onsager-Feynman quantum vortices 
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has been modeled by Feynman 1 and others (compare study by Ladner 2 ). For instance, at the exit 
of a tube the shedding of "smoke ring"-like vortices has been envisioned as the critical process. 
The resulting critical velocity is proportional to the quantum of circulation (h/m = 10 -3 cm 2 /s), 
and it is inversely proportional to the tube radius (D/2). Very small critical superfluid velocities 
result for radii larger than 1 cm. 

Another critical flow rate is predicted by the Piteavskii-Mamaladse 3 version of the Lan- 
dau theory of higher order phase transitions 4 . Starting from basic superfluid parameters, it has 
been possible to quantify critical velocities of superfluid orifice flow 5 


Critical Counterflow. A different class of criteria pertains to counterflow of superfluid and nor- 
mal fluid, involving a finite heat flow rate, at zero net mass flow (ZNMF). A special case has 
been the Dimotakis criterion 6 relying on the mutual friction parameter of the Gorter-Mellink 
transport concept for counterflow 2 . The present simplified view of this criterion ignores the vari- 
ous transition steps known for ducts with adiabatic walls (ZNMF ducts). Details of the transition 
have been studied and surveyed 7 ’ 2 . In the present porous media however, a finite pore size distri- 
bution masks the "fine structure" of the transition. Therefore, the present simplified approach 
ignores these details concentrating instead on steady values (usually obtained for externally 
imposed decreasing heat flux densities (dq < 0)). In this manner, hysteresis effects, leading to 
large critical rates at vortex shedding onset, are de-emphasized. The Dimotakis result for the 
critical counterflow velocity ((w c = | v c -v s | c ) is written as 


w c = (l/7C)/[DpA G M,c] 


( 1 ) 


(D duct diameter, p liquid He II density). The mutual friction between vortex lines generated 
from superfluid and normal fluid excitations is described by the Gorter-Mellink function A GM , a 
function of temperature T and pressure P. The critical value A G m, c has been documented only 
slightly. However for the asymptotic regime of chaotic counterflow, with size-independent q- 
values, rules of similitude have been introduced 8 . The Gorter-Mellink function of this similitude 
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approach for the high-q range is formulated in terms of basic thermophysical properties of He II 


Agm = Kgm(P/Ps)/ti„ 

(2) 

(Kgm = const = 1 1.3, T| n shear viscosity of normal fluid). Starting from this high-q information, 
an order of magnitude estimate of the critical counterflow velocity (w c ) is considered: the high-q 
function for Agm is substituted for the critical value at low q = q c . One then arrives at 

w c = (K^ M M)Tln(Ps/P)/[Dp] 

(3) 

In the discussion of equation (3), it is noted that at low T we have a very small normal fluid den- 
sity [p n < p; (p s ~ p)]. Furthermore, the critical counterflow speed approaches the critical nor- 
mal fluid velocity (w c = v,*). Several authors have addressed this condition, noting that a normal 
fluid Reynolds number criterion ought to become valid as soon as the normal fluid speed 
increases significantly. According to the Reynolds number criterion, the total density has been 
used (p = ps+pn) with reference to normal fluid values in speed and shear viscosity 

Re c = v„ c pD/n n ~ w c pD/n n 


Pn<P 


(4) 


The corresponding counterflow criterion is then recast, paralleling the normal fluid Reynolds 
number. From equation (2) we obtain 


Re c ~ (KaM/nXPs/p) 


(5) 


The constant (K^mM) turns out to be 459. In other words, this constant approaches the same 
order of magnitude as the classical Reynolds number for Newtonian fluid flow through ducts 
(Re c = 2300). Thus, the simplest estimate, permitted by the present approach, appears to be a 
combination of statement (4) with the temperature dependence of equation (5). The resulting 
Reynolds number for normal fluid flow, using the total density, is written as 
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Re c = Re^ps/p) 


( 6 ) 

The heat flux density at zero net mass flow results from transport of (pST) per unit volume at the 
speed v n 


q = pSTv n 


(7) 


(S entropy per unit mass of liquid). Thus the normal fluid Reynolds number (4), at critical con- 
ditions, is equal to q c D/(r| n ST). Figure 1 presents data which provide only limited support for 
equation (6). 


Critical temperature gradient condition. The dimensionless form of the critical T-gradient is for- 
mulated assuming a porous media geometry. The "classical" flow of Newtonian fluid at low 
speed has been described by Darcy’s law. The mean speed of the fluid is conveniently expressed 
as approach velocity v 0 (= superficial speed). Zero net mass flow studies at low temperature 
have led to the adoption of an analog of Darcy’s law for the normal fluid. This particular equa- 
tion has become quite useful for the characterization of several flow rate conditions 9 . The Darcy 
law analog is written in terms of the approach velocity for normal fluid flow (and heat flow). 


v n o = K„|gradPrl/nn 


( 8 ) 


The driving gradient is the thermomechanical pressure gradient, described by London’s equation 


grad P = p S grad T 


(9) 


The normal fluid permeability is K n . The Darcy permeability of a clean porous plug, without 
impurity condensate layers, ought to be very close (within 1%) to the room temperature value for 
Newtonian fluid, unless stress concentrations and other effects change geometric conditions. 
After insertion of equation (9) into (10) one may obtain a dimensionless form of Darcy’s law 
analog, including the dimensionless T-gradient. Both sides of equation (8) are multiplied by 
pLcMn- Lc is the characteristic length of the porous medium (with L } = K n ). Because of the 
correspondence in terms of permeability, L c is a measure of the fluid throughput capability of the 
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porous medium. The dimensionless result obtained is 


pS | grad T | L^pAln — v c pLc/r| n 


( 10 ) 


The left hand side of equation (10) represents the dimensionless T-gradient, and the right hand 
side is again the normal fluid Reynolds number. At critical conditions the subscript c is intro- 
duced. One possible postulate for a quantification of critical counterflow is the condition of a 
superfluid density-dependent function. This implies a temperature-dependence with the con- 
straint that Re c approaches zero, as T approaches the lambda point (Tx)» while (Ps/p) decreases 
toward zero. 


Another postulate, supported by literature data, is the adoption of a temperature- 
dependent (or order parameter-dependent) function of the dimensionless critical T-gradient, 
called Nvt c - 

NvT c = pS|gradT c |pL c W 

Investigations of porous media 10 and of a duct 11 suggest that Nvt c varies from zero, at Tx. to 
finite values, proportional to (Ps/p), when T is lowered. Thus, there exists a critical number (N c ) 
such that the simplest counterflow criterion for a specified porous medium is 

Nvr e (Pn/ps) = N c 

(12: 

The value of N c is to be determined experimentally for a particular porous medium as long as 
further porous plug parameter effects are unknown. 
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Experiments 


Experimental porous media studies have been conducted with sintered stainless steel 
plugs. The investigations have been quite similar to ZNMF studies of other porous media, e.g., 
reference 12. A chamber equipped with heater and thermometer is separated from a main He II 
bath by means of a porous plug (inset of Fig. 2). Other chamber walls are insulated. The wall 
material is phenolic (outer diameter of Fig. 2 [inset] is equal to 2.5 cm, overall length is 3.6 cm). 
The stainless steel plugs employed have the following dimensions: PKS plug, very thin disk 
with a diameter of D = 0.64 cm, thickness L = 0.075 cm; PSM plug, cut from a large disk, diam- 
eter D = 0.79 cm, L = 0.175 cm; PSM-type plug of small dimension, D = 0.323 cm, L = 0.16 cm. 
The latter plug has been arranged as a triple set, acting as three "flow resistances" in series. The 
sintered stainless steel plugs are mounted with axis oriented vertically. The heater inside the 
chamber has been wound using constantan wire (75 |xm diameter) to approximately 100 ohms. 
The carbon resistance thermometers inside the chamber produce a thermogram, with non-linear 
T-scale recording the thermometer voltage versus time t. 

The liquid He II inside the chamber can communicate with the outer liquid reservoir, per- 
mitting a very small mass inflow during initial heating. Heat removal by normal fluid, i.e. 
outflow out of the chamber, takes place during the ZNMF experiment. In contrast to previous 
studies, 12 a modification has been introduced. A transient technique is employed imposing an 
externally applied step power input to the heater inside the chamber. 

Figure 2 is an example of the thermometer voltage versus time record at relatively low 
heater power. At moderate and high power inputs the thermograms obtained show additional 
features, such as onset of a phase transition, oscillatory modes of the temperature (T) versus time 
(t), and possibly large T-excursions into the vapor regime. The various phase transitions have 
helped in the diagnostics of the records, in particular identification of the critical rates. The final 
T-differences for steady transport, seen on the thermogram at long times, have been the basis for 
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the determination of critical parameters. At each bath temperature, the critical parameter set 
(AT C [Q C ]) separates the linear transport regime (= laminar normal fluid flow) from the non-linear 
regime. 


Small power inputs are necessary to cover the linear regime of (absolute) heat input Q 
versus temperature difference. Above the critical value, a non-linear Q(AT)-function is observed 
(e.g. Fig. 3). Figure 3 plots the approach heat flux density Q/A = q 0 (A total plug cross sectional 
area). For most plugs, one obtains readily the permeability of normal fluid flow (k^) from the 
records, at least at low T. The defining equation (8), in conjunction with equation (9), may be 
expressed in terms of a property-dependent function F(T) = [pS/rjnKpST). We make use of the 
ZNMF equation for the approach heat flux density. From q = pSTv no , and Q = Aq 0 , one obtains 


Q - A(pST)K T1 pS[AT/L]/'n n 

After insertion of F(T), the normal fluid permeability is expressed as 


(13) 


K n = [Q/AT]LA/F(T) 


(14) 


At low T, the function F has very small values rising rapidly as T is raised. Therefore, for a 
specified heat input Q, the AT-value in the linear range may be relatively large. In contrast, at 
"high" T, near the lambda point, only a small AT is established for the Q- value specified. The 


inset of Figure 3 shows F(T). 


Discussion of Results 


Critical temperature difference function. Plug data obtained for long times as steady value 
AT(Q), provides a critical parameter set at the intersection of the linear characteristics with the 
non-linear function. Figure 4 presents the quantity Nvt/Lc as a function of liquid bath tempera- 
ture. It is noted from equation (11) that the normalized critical T-difference involves the 
geometry only via the plug thickness (L). Further, there are only thermophysical He II proper- 
ties aside from the AT c -value. In Figure 4 all results obtained, within data scatter, and for the 
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temperature range covered, appear to be consistent with a variation of [NvTc/Lc ] proportional to 
the superfluid to normal fluid ratio (ps/p n )- It is noted that the density ratio has the value unity 
near 1.95 K. Therefore, an appropriate constant [NvT c L(T 3 Pn/Ps] may be read off directly from 
the log-log plot in Figure 4. The function [10 11 (Ps/Pn)] in cm -3 , has been included in Figure 4. 
Previous data 11 have been found to be proportional to (Ps/Pn) during early studies of a porous 
medium of high Darcy permeability. 

Critical Reynolds number. The critical Reynolds number concept has been discussed primarily 
in early He II work using ducts. 2 In Figure 1 the critical Reynolds number equation (6) has been 
employed for a comparison. A reference critical Reynolds number of 2300 has been inserted. 
Various early data sets show little agreement with each other, presumably because of different 
techniques employed and because of different observation limits for critical conditions. How- 
ever there are common tendencies. As T is lowered, starting from Tx, there is a monotonic 
increase of the critical Reynolds number with (Tx-T). At low T, all of the data appear to indi- 
cate a tendency to increase less drastically in the (absolute) derivative | dReg/dT | . In fact, some 
data sets appear to level off at Reynolds numbers on the order of magnitude 10 3 . Results by 
Chase 13 document this behavior. 


Figure 5 is a plot of data 13 for four tubes with diameters which range from 0.08 cm to 0.4 
cm. The inspection of Re c provides support for the T-gradient criterion near the lambda transi- 
tion. At low T however, the critical Reynolds number criterion appears to be a constraint. Thus, 
at this time we propose a simple empirical interpolation in the following form 


Re = [Re r ^ + (p n /{p s 550))- 4 ]- 1/ « 


( 15 ) 


This function is included in Figure 5. It is seen that equation (15) describes the critical Re c 13 - 
data quite well. 


Conclusions 


32 



It is concluded that critical counterflow in porous media has a temperature dependence 
supportive of a dimensionless critical temperature gradient This quantity is found to be propor- 
tional to the superfluid to normal fluid density ratio, for the data sets obtained in the present 
work, as far as the T-dependence is concerned. 
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Figure 1 Critical Reynold's number versus T; (a), function Re € » 
550 (pjp^ (b) # Equation (6). O. Brewer-Edwards. 108 m; +. 
Chase, d -0.159 cm; A* Griffiths «f a/.; A A. Vinea 0.35 cm. 0.53 
cm; ■, Staas et *L 0.255 m 



Figure 2 Thermogram: temperature T versus time; plug PA/FA 
PSM OIS-tr-3.2 x 1.6. Inset: schematic diagram of experimental 
apparatus with triple set of sintered stainless steel plugs. H. heater; 
PH. phenolic sleeve; PP. sintered Porous plug: TH. carbon ther- 
mometer 


36 







ORIGINAL PAGE IS 
OF POOR QUALITY 



*r oo 


Figure 3 Approach hast flux dansity q 0 versus T at Ta*™- 
1.91 K for plug K10S02-6.4 x 0.75. Permeability. x - 2.06 x 10~'° 
cm 2 . Inset: F( 7) of Ho II properties; Equation (14). 
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Fig. 4. Critical quantity N^j,^/L c 3 an a function of temperature; 

A plug K 10 S 02-6.4x0.75 ; (2.2) (p 8 /p n ) i 0 13 

• plug DO 01-S-05-7.9xl.75 ; 7.5 10 12 (p g /p n )’; 

o plug PA/FA PSM 01 S-tr-3.2 x 1.6 ; 1.3 10 12 (p s /p ) ; 

in * rn 

T porous medium , Ref. 10 ; 2.2 10 • 
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CHARACTERISTIC PARAMETERS RELATED TO SUPERCONDUCTOR-COOLANT INTERACTION 
T. H. K. Frederking 

University of California, Los Angeles, California 90024-1592, USA 
ABSTRACT 

In the area of basic heat transfer mechanisms of Helium heat transfer and related influence on stability of superconduct- 
ing magnets, thermal boundary conditions are among the sets of influential constraints. Important parameters, in particu- 
lar characteristic lengths, are considered for the two sub-systems of a composite, i.e. solid state domains and fluid 
domains, noting examples of high transition temperatures. 


Paper presented at the JOINT SEMINAR (US-Japan) on 
BASIC MECHANISMS OF HELIUM HEAT TRANSFER AND RELATED 
INFLUENCE ON STABILITY OF SUPERCONDUCTING MAGNETS, 

Aug. 29 - Sept. 2, 1988 , Fukuoka , Japan, paper of Session IB 
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INTRODUCTION 


The evolution of large superconducting magnets during the last 25 years has been based on the composite NbTi-Cu, a 
most popular metal superconductor (MSQ. Its transition temperature T c of about 9 K requires Helium-4 as coolant 
The advent of high T (ceramic) superconductors (HTSC) in the last two years has posed questions concerning stabiliza- 
tion approaches (e.g. 1 ) developed for Helium coolant systems. It is the purpose of the present investigation to inspect the 
conditions of superconductor-composite stability. The emphasis is on thermodynamic and transport phenomena related 
to the sizing of stabilizers. First solid state equilibrium is considered, subsequently fluid domains for both Newtonian 
and superfluid (non-Newtonian fluid) are inspected. Finally examples of interactions between composite and coolant, 
involving solid and fluid parameters, are presented. 

CHARACTERISTIC SUPERCONDUCTOR PARAMETERS 


The thermodynamics of a superconducting system is referred to using simplified assumptions associated with "supercon- 
ducting electron fluids” and with "normal electron fluids". The extended Ginzburg-Landau (GL) (e.g. 2 ) thermodynamic 
theory is the basis, along with Tuyn’s law, for the critical magnetic field, written as flux line density 


B c = B„(l-t 2 ) 

B 0 zero Kelvin value; t = T/T c ; T temperature. The related "two-fluid" number density is 


( 1 ) 


n, = n 0 (l-t 4 ) 


( 2 ) 


n, = number density of electrons participating in the ordering process leading to low lying superconducting energy states. 


There are three important lengths: 1) the electronic mean free path (MFP) e ; 2) the coherence distance (£); and 3) the 
(weak field) penetration depth (Ap) for the exclusion of magnetic flux from bulk material (Meissner effect). There are 
two basic types of matter 3 : type I superconductors characterized by a ratio k = (Xp/£) below unity, and type II supercon- 
ductors with large k values. Type I is characterized by the thermodynamic critical field (B c of Eq. (1)). Type II material 
is characterized by two critical magnetic fields (B c i and B C 2 ). At the lower critical magnetic field (B c i), quantum flux 
lines start to penetrate into the superconducting bulk material. At the upper critical magnetic field (B C 2 ) flux line thread- 
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ing of the bulk material is complete such that the normal state is restored. MSC materials used for magnets have k 
values on the order of magnitude of 100. The £ values approach the order of the relevant lattice constant, i.e. values of 
atom spacing in the lattice. HTSC appear to show similar features with additional constraints of anisotropy. Neverthe- 
less, it appears to be possible to arrive at estimates also for these systems assuming significance of the extended GL 
theory (Appendices A and B) for magneto- thermodynamic order of magnitude estimates. The parameter set comprises 
B c and the number density (n,) of the condensation process. 


Figure 1 displays the critical current density at zero field expressed as ratio (jc/B 0 ), for a transition temperature of 90 K. 
The thermodynamic result for B C 2 = V2 kB c is (Appendices B and C) 4 - 5 (n« = number density of T = 0 K) 


jo = [Bc 2 /(V 2 K)][2/3]^2e[n 0 /(m e Mo)] V4 

(3) 

With MSC materials, high jc values have been realized with suitable thin films. For HTSC it appears to be reasonable to 
expect potentially useful results for similar films of high quality with flux pinning. The anisotropy of the HTSC materi- 
als however introduces different number densities (n,) in different directions, causing anisotropic critical current densi- 
ties j c . 


The simplified ideal thermodynamics leaves open considerable questions about stability of flux lines and stability with 
respect to dissipative events elsewhere (e.g. boundaries of coolant domains subject to phase transitions). There is field 
penetration into magnet material during magnet energization. During fast flux line motion, adiabatic conditions may 
prevail. Interdependence between critical current and line motion requires rate control by suitable pinning conditions. 
Because of normal states in vortex cores, and because of the necessity to bring in flux lines into the superconductor, dis- 
sipation is present. This causes magneto thermal heating influenced by temperature-dependent effects: the critical 
current density decreases with T while the heat capacity increases. There is the potential for adiabatic flux jumps 6 with 
final T > T c . Two major steps in the magnet developments have been stabilization using normal metal, ("N" = Cu, Al) 
as "shunts”, and subdivision of the superconductor into fine filaments 1 . A simple energy criterion for "magnetic pres- 
sure” thermalization leads to an upper bound for the filament diameter (dF> (Appendix D) 


d F = { PC Vjc)Po 1 (2MopC p AT]' /i 


(4) 
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For HTSC 6 , with large enthalpy, cfc becomes rather large for moderate current densities. Large enthalpy differences, ZpAT 
= [H(T c >-H(Tb)], appear to be available for hydrogen-cooled HTSC 7 , if safety concerns are secondary (p density). 


In addition to filament concerns, the MSC magnet R&D work has relied on well defined cooling systems. The related 
heat transfer has been considered early in conjunction with normal metal shunts introduced by Laverick and Stekly 1 . 

The question for optimal cross sections has stimulated modeling of the scenario of a local quench. Relying on a stability 
parameter (a), the Joule dissipation of steady operation (''cryostadcs") has to be proportional to the heat removal capa- 
bility of the system. Fourier diffusion ought to occur from T c down to the vicinity of the bath temperature (Tb). 
Governing parameters are the heat transfer coefficient (h), the stabilizer cross section (Ac), the fluid-wetted circumfer- 
ence (C c ) and the electrical resistivity of the stabilizer (p e ). The limit of stable operation is reached for a on the order of 
unity. The required area is 


Ac^PcLVlhCcCTc-TB)] 


( 5 ) 


(I c = critical current). If h is large, excessive cross sections can be avoided. If the fluid undergoes unfavorable phase 
transitions, leading to a low h, rather large Ac, and along with it low overall electrical current densities are reached. This 
enlarges the volume of a magnet specified for a set of constraints. 


FLUID DOMAINS: NEWTONIAN FLUIDS 


The abundant literature on fluid convection permits numerous order of magnitude assessments for steady transport. The 
heat transfer coefficient (thermal conductance h) is proportional to the Nusselt number (Nu). The latter is proportional to 
the thermal conductivity (k) which in turn is proportional to the specific heat. Other parameters for a specified flow 
regime may exert secondary effects (possibly large in terms of pressure drops). As pointed out by Reynolds, pressure 
drop and heat removal rates are coupled, presenting challenges to coolant duct optimization. 

Forced convection. The characterization lengths for short winding elements may be laminar entrance lengths, e.g. 
boundary layer thicknesses (8,8 t) whose ratio depends on the Prandtl number (Pr = qcp/k), q shear viscosity. Hydro- 
dynamic boundary layer thickness 8 


43 


6 = 5[tix/(vp)] 


( 6 ) 


Thermal boundary layer thickness 


§r = 5/NPr 

(7) 

(v velocity, x distance from entrance of rather wide coolant duct). In turbulent fluid, vortex shedding involves large bun- 
dles of quantum vortices whose discrete nature is of no concern in "eddy distributions". The steady state vortex distribu- 
tion of well developed turbulence has been characterized in terms of layer models (e.g. wall domain with laminar flow 
features, buffer layer). 


Natural convection modes. The natural ("free convection") mode may appear to be attractive at first glance in part 
because of die lack of a pump requirement. However, h on the order of 0.01 W/(cm 2 K) is rather poor causing great 
concern upon inspection of Eq. (5). Thus, the loss-of-flow scenario, starting from a faced flow system, is an unpleasant 
one. Nevertheless, the possibility of high h values near the thermodynamic critical point (P c ,T c t) has prompted research 
of phase transitions without latent heat above and near P e . Because of He I and He II, there are numerous phase transi- 
tions in Helium-4. For HTSC in nitrogen, the supercritical states of forced flow are quite a bit away from the transposed 
critical curve (TCC). Thus, numerous phase transition possibilities of He near T c of the MSC appear to be of no serious 
concern for HTSC. 


An exceptional case has been the centrifugal buoyancy-induced free convection, e.g. in a SC generator for elec- 
tromechanical energy conversion. Coolant duct boundary layers have a thickness 


with an "effective pressure gradient" 


St - [(Lgkq)/(pc p | grad P g 1 ] v * 


( 8 ) 


I grad P g | = |a(p~-p)| 


(9) 


L g is the length of winding package exposed to acceleration (a), e.g. ”g”, or centrifugal buoyancy rQ 2 (r radius, Cl angu- 
lar frequency). For power station size a.c. synchronous generator windings, the superconductor may be cooled con- 
veniently using the available gradient with gravitational acceleration (g) replaced by tCI 2 . The order of 1000 g may be 
reached causing thinning of the boundary layers, Eq. (8). Coefficients h comparable to nucleate boiling coefficients may 
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be reached. 


Phase transitions. The most popular case has been nucleate boiling with h on the order of 1 W/(cm 2 K) in He I. In con- 
trast, the high T excursion leading to film boiling, is characterized for liquefied gases in pool situations by h values on 
the order of 10* 2 W/(cm 2 K). The drastic lowering of h, by a factor of 100, has caused sufficient stability concern. 
Therefore, a different scenario has been formulated for the nucleate pool boiling cryostatics. The parameter (a0 
assumes the dual role of a safety factor and a stability number. For the maximum heat removal rate at the "peak" heat 
flux density (qp) one ought to reach the highest current density if cti reaches unity. The stabilizer cross section is 


Ac i cti 'IcVtqpCc) 

( 10 ) 

There is a more conservative approach using the "recovery" value q„u„ at the minimum heat flux density of the boiling 
curve. The simplification of the pool boiling curve suggests replacement of qp in Eq. (8) by q„un (and possibly replace- 
ment of cti by another value 012 ). A rather drastic enlargement of the area A« results from this procedure. A middle-of- 
the-road approach appears to be offered by the "equal area criterion" 1 , adopting a mean value of the boiling curve heat 
removal rate with respect to the transition temperature T c . 


FLUID DOMAINS: SUPERFLUID He II 


Bulk He II has a low thermal resistance and excellent transient heat removal capability. Therefore, surface domains, 
such as insulating layers, and the Kapitza resistance are "seen" easier than in He I. London has pointed out the convec- 
tive nature of He II heat transport while recent research emphasizes vortex system behavior 9 * 10 . Quantum vortex bundles 
produce "mutual friction" retarding the flow of "normal fluid" which, in two-fluid model language, carries the entropy 
and the thermal energy of He I. Superfluid produces the thermomechanical effect, a pressure increase useful for ther- 
mally induced "forced flow". There is the benefit of independence from gravity force orientation. 

Sizes of vortex bundles ought to affect the flow in narrow cooling passages. Phase separator data for space systems, 
e.g. 11 and zero net mass flow modes in similar porous media have shown that passages below the order of magnitude 
10 -3 cm do not permit high transport rates, predicted for bulk liquid, and desired for large He II-cooled magnets. This 
size influence has been represented by a fluid "rate constant" (Kvlps, Kznmf). The latter is the ratio of an effective 
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counterflow speed, v e = [q/(p,ST)] divided by the speed predicted from mutual friction of bulk He II (S entropy per unit 
mass). 


K* = v e /[S|VT|(Vpn)(p^p )] 1/3 

Figure 2 presents rate constants versus the Darcy permeability k of a porous medium. 


( 11 ) 


There appears to be an extraordinary richness in phase transitions at the low temperatures of He II. This may produce 
peculiar loss-of-flow quench scenarios. There are very attractive features related to the possibility of consistently high h 
values, a large enthalpy difference with MSC, a "heat pipe" like transport rate in the liquid for not too narrow cooling 
passages. In addition, the thermomechanical effect may be used for pumping of liquid through magnet windings without 
excessive mechanical (electrical) power taken from an external device. 

SUPERCONDUCTOR-COOLANT INTERACTION 


As illustrated by various stability parameters, there appear to be numerous possible approaches referring to the thermal 
boundary conditions in the form of h. Often, simplified boiling curves for pool geometries have been adopted in order to 
assess the "degree of stability". Experimentally, exceeding the limits ought to lead to a quench with a finite propagation 
velocity. In this regard not only the criticism of unrealistic heat removal functions but also the danger of arbitrary stabil- 
izer area (A) enlargement appears to be worth noting. Minimum propagation zones (MPZ) (e.g. Ref. 12 of this seminar), 
may function as onset of quench spreading. However, for film boiling phase transitions, a stationary "hot gas" domain is 
possible while dissipated energy is conducted away by Fourier diffusion. This detrimental case is characterized by a 
thermal decay length of Ld = [kA/(hC c )] H - (Dk/h) % of relative large extent (D diameter of composite conductor). A 
possible situation of the MPZ is a minor T increase within the quenched region, a relatively short current sharing domain 
and a major decay over the length Ld on both sides of the quenched domain. 

Unsteady energy deposition criteria permitted by the coolant appear to bring less conservatism and more attractive 
overall current densities to a magnet geometry. In conjunction with a.c. conditions, the lack of a universal coverage of 
superconducting equipment stability has been noted (Ref. 13). 
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In the present context a simplified example is considered using "nucleate cooling stability 14 - 15 for HTSC cooled by He I, 
liquid hydrogen and liquid nitrogen in nucleate boiling modes. The peak value qp is the characteristic parameter for the 
coolant limit. The electrical resistivity (p e ) of Matthiessen’s rule for Cu is the dissipative term’s characteristic property. 
A residual resistivity ratio of SO is assumed in conjunction with a simplified concentric cylinder conductor geometry (r = 
D/2; 2r e = superconductor diameter). The overall current density limit is obtained from aj = 1 as 


T= [(2q P /p e )(r 2 -r e 2 )/ r 3 ] ' A 

Optimization of j with respect to r, for a constant r c , results in a maximum of j =J* 


( 12 ) 


? = [4qp/(3 Pe r)]* 

(13) 

Figure 3 is a plot of J* versus T using peak heat flux densities of nucleate boiling curves with emphasis chi Kutateladse- 
Zuber functions 16 . It is seen that there is a lack of a drastic order of magnitude change as one switches from one coolant 
to another one potentially useful for HTSC. The J* values are on the order of magnitude 10 4 A/cm 2 near maxima of the 
J*(T) functions. This is in contrast to the potentially large intrinsic critical current densities, in zero field, of the thermo- 
dynamic model equation (3). In the presence of premature flux jumps, models for these phenomena 17 may be useful. 
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Fig. 1. Intrinsic critical current density expressed as ratio jc/B„ versus temperature ratio 
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Fig. 2. Mutual friction "transport rate constants" K/ (j = VLPS, ZNMF) as a function of the Darcy permeability for 
porous media; 

^ zero net mass flow mode 

® VLPS vapor liquid phase separation mode. 
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Fig. 3. Transport parameters of nucleate boiling and superconductor/boiling liquid interaction: 

a. Peak heat flux density of nucleate boiling of saturated liquid; 

b. Optimum overall critical current density of composite; radius D/2 = 0.2 cm; RRR * residual resistivity 
ratio 50; Cu. 
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APPENDIX A. MAGNETIC FIELD PENETRATION DEPTH (WEAK FIELD CONDITIONS) 


Consider a boundary between "vacuum" and a type I superconducting state. The applied magnetic field is parallel to the 
surface of the superconductor. The shielding currents associated with the Meissner effect decay in a domain of finite 
thickness, and so does the field inside the superconductor. The energetics of this situation is considered minimizing, in a 
simple model, the total energy, with respect to the distance L from the surface at z=0. One contribution is the magnetic 
pressure-related energy, per unit area, of AP m L = LB 2 /(2(Xo). The other important contribution is the electron kinetic 
energy of m e n,v 2 /2; n, number density; m* electron mass, v velocity; where j* = n,ev is the electrical current density. 

The latter is related to B, applied at z=0. From curl B = |i<j e we obtain the dominant contribution: j i ~ [B 2 /(po)]/L. 
Minimizing the sum of the two contributions, dEuJdL -> 0, one obtains the magnetic field penetration depth 

X p =[m«/(n g e 2 Mo)] ,/4 


(A.1) 

It is seen that the present energetic assessment contains only one parameter, e.g. n, of Eq. (2) which influences the pene- 
tration depth. 
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APPENDIX B: SIMPLIFIED MAGNETO-THERMODYNAMIC SUPERCONDUCTOR MODEL 


Consider a macroscopic superconducting system. The Ginsburg-Landau (GL) phase transition equations (e.g. Ref. 2,4) 
are used and extended to low temperatures. Phenomenological equations for the macro-wave function, with number 
density | y | 2 ~ n,, have equilibrium state minima with respect to normal states. The free energy density is g = 
g n -a I V | ^ (P /2) | y | 4 . The equilibrium state free energy difference, between normal and superconducting states, is 
obtained as 


Ag = n,f(t)oto 

(B.l) 

Using similarity functions, e.g. Tuyn’s law for the critical magnetic flux density, one arrives at simple functions of t 
(t = T/T c ). At T = 0 K, the GL functional has values 


cto = 2APoAio 


(B.2) 


and 


Po= Oto/llo 

(B.3) 

(AP 0 = B^/(2po) = magnetic pressure difference, equal to the maximum free energy density difference). The number 
density no = |y<>ol 2 is quite low for commercial metal superconductors as only the fraction (Tc/Tp) participates in the 
ordering process ("condensation" of the Bose-Einstein kind in momentum space) [Tp Fermi temperature]. The tempera- 
ture effects on the functions (tx/Oo) and (P/p 0 ) are shown in Fig. B.l. Figure B.2 illustrates the lowering of the free 
energy density, expressed as Ag/Ago versus [ | y 1 2 ] Vi . 


Concerning characteristic lengths, the following functions are obtained: the penetration depth (Appendix A) varies as 


y**= tMt 4 )]-* 


Noting the GL definition of the coherence length £ =fi/[2m e a] V4 , we have 


(B.4) 


^o=[(l+t 2 )/(l-t 2 )]‘ A 


(B.5) 


The term (a/m e ) may be expressed as kinetic energy. Considering <£v-dl = NT 0 around a quantum flux line (vortex line), 
we make use of the quantum of circulation T 0 = (2rth/m) (!i Planck’s constant divided by 2n). The energetically favored 
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lowest quantum state, with number N=1 is inserted in the energy per mass: 


(o/m e )=-F 0 2/(2^2) = v2/2 

(F 0 = TJ2n) flux quantum <|> 0 = 2rfli/(2e) = 2.07xl0" 15 Wb = 2.07xl0- 7 G cm 2 . 

Physical constants 

e = 1.602 10- 19 As 
me = 9.1095 10~ 31 Kg 
|Xb = 4jc 10* 7 £2s/m 
11= 1.054 10 -34 J-s 

Reference value: iv> = 1 m -3 : 

Oo/Bo) = 0. 1630 (A/m^/T for T = 0 K 


(B.6) 
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Fig. B.2. Free energy density difference ratio versus order parameter ratio; temperature ratios T/T e : a = 0.9; b = 0.8 
c = 0.6; d = 0.4; e = 0.2; f = OK. 
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APPENDIX C: SOLUTION FOR "UNIFORM FLOW' 


The "intrinsic" solution for uniform flow (Eq. 3) derives from inclusion of the kinetic energy (me/2)v 2 in the free energy 
functional (v velocity of the "superflow” of electrons). 


g = g-a | \|t | 2 + (P/2) | v | 4 + (m«/2)v 2 

(C.l) 

The free energy is minimized with respect to | y | 2 for a constant temperature system at uniform flow. Further 
j e - (2e) | y 0 1 2 v has a maximum value of j c . Aside from the number density and constants, only thermodynamic quanti- 
ties appear. Thus, an "intrinsic" value is reached. The difference Ag versus v has a point of inflection, leading to a situa- 
tion similar to "choked nozzle flow". The limiting rate, Eq. (3), may be written in terms of T 0 as 


jc = (4/3)3-*e|y 0 |?Fo£ 

(C.2) 

The energy per mass (a/m e ), Eq. (B.2) of Appendix B, relates to the thermodynamic magnetic pressure difference, i.e. to 
B c . For T=0 one obtains 


(ct/m.) =-F 2 /(2Jjo) = 2AP 0 /(Voo) 2 m. = (B ^/^/((Voo) 2 ™.) 


Eliminating T 0 from Eqs. (C.2) and (C.3) one arrives at 


jo — Bo 


2 

7 


3/2 


2e 


no 


m«Jio 


vi 

2 * 


The temperature dependence of the critical current density is obtained readily (Appendix B). 


jc/jo = (l-t 2 )[l-t 4 ]’ A 


(C.3) 


(C.4) 


(C.5) 
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APPENDIX D: FILAMENT SIZE GUIDE LINES BASED ON SUPERCONDUCTOR ENERGETICS 


Consider a thin slab ({filament model") of superconductor of width dF with sufficient pinning centers. An 
externally applied field B produces an energy (per unit volume) proportional to the magnetic pressure difference 


AP m = B 2 /(2po) 

(D.l) 

Once a limiting condition with dissipation on the order of pp c (per unit volume) is attained, the enthalpy difference AH 
has to be sufficiently large to prevent "take-off" to high normal state temperatures: 


J^e-ptHCTehHCTo)] 

(D.2) 

where T c is the transition temperature and T 0 the superconductor temperature, assumed to be close to the coolant tem- 
perature. The enthalpy per volume pAH may be expressed in terms of the specific heat (per unit volume, of pep), and in 
terms of the temperature difference T e -T 0 = AT. 


•A «A 

According to Maxwell’s equation, curl B = p<j, there is a current density on the order of 

j - 

(D.3) 

(D.3) is related to pinning models (e.g. Bean model) when j * je is the critical current density of the limiting pinning 
strength. Equating the energy densities (D.l) and (D.2) one arrives at a field value (order of magnitude) 


B - ((2|i 0 )pC p AT] v * 

Making use of (D.3), for the critical value, and of (D.4) we obtain a simplified limit for the filament of 


(D.4) 


dF = p(l/jc)Mo-U2McP^T]» 

(D.5) 

The factor f) is on the order of magnitude unity. A metal superconductor is characterized by the following orders of 
magnitude: 

Specific heat per unit volume: pep - 10 3 J/m 3 K) 

Temperature: AT - 1 K 
External field value: B - 1 T 
Critical current density: j c - 10 9 A/m 2 

The resulting filament size turns out to have the order of magnitude 10 pm. 
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APPENDIX 


D 

ABSTRACT OF THE THESIS 

The Characterization of Sintered Stainless 
Steel Porous Plugs for the Zero Net Mass Flow 
Mode of Helium II Including Nonlinear Flow 

by 

Parimal K. Khandhar 

Master of Science in Chemical Engineering 
University of California, Los Angeles, 1989 
Professor Traugott H. K. Frederking, Chair 

The Zero Net Mass Flow (ZNMF) transport mode of He II 
has been investigated. In this mode the normal fluid flow 
emitted from a heater is compensated by an inflow of the 
superfluid into the heating chamber using two-fluid con- 
cepts. Three sintered stainless steel plugs employed have 
Darcy permeabilities on the order of 10 1 ° cm 2 ; (10 10 ae: 
4.44 +_ 0 . 03 #1, 2.06 +, 0.15 #2, 9.22 +. 1.0 #3). Because of 
limited data in the nonlinear ZNMF regime, heat flux den- 
sities (q) have been obtained as a function of temperature 
difference (AT) in this regime. The heat flux density is 
found to be permeability-dependent, f(ae) = K z N „ F = 

[ ( 1/Ko „ ) 2 + (ae r „ £ /ae) ]" 1 ' 2 , similar to the vapor liquid 
phase separation results of Yuan. The following semi- 
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empirical equation is proposed to fit nonlinear results. 
q= [ ( l/K, n ) 2 + (« ro£ /<s)]~ 1/2 f S 4 T 1 i lf IVT I /p ] 1 ' 3 

ae r . £ = 3 x 10' 10 cm 2 ; AT << T 

(S entropy, density, p j (j = n, s : subscript n = normal 
fluid, s = superfluid), = shear viscosity, K r ,„ = Gorter 
Mellink mutual friction constant = 11.3 for wide ducts. 

Two plugs tested are found to have K* N „ F = 0.3 +, 0.01 (#2) 
and 0.556 + 0.04 (#3). An alternate equation has been de- 
duced removing a conceptual difficulty of the above semi- 
empirical equation. 

q = K 4 /(L L c ) 1 ' 4 [fj S 5 T 4 |VT|/ 9 * 9 „ ] 1 ' 4 

(L = thickness of a plug, K 4 = Heat transport rate con- 
stant; 3.1 +, 0.2 (plug #2). 

The present work has been extended up to the vapor- 
ization limit. The data exhibit superheating of He II 
across the vapor pressure curve. 
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